Experimentally witnessing the quantumness of correlations 
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The quantification of quantum correlations (other than entanglement) usually entails laboured 
numerical optimization procedures also demanding quantum state tomographic methods. Thus it 
is interesting to have a laboratory friendly witness for the nature of correlations. In this Letter 
we report a direct experimental implementation of such a witness in a room temperature nuclear 
magnetic resonance system. In our experiment the nature of correlations is revealed by performing 
only few local magnetization measurements. We also compare the witness results with those for the 
symmetric quantum discord and we obtained a fairly good agreement. 

PACS numbers: 03.67.-a, 03.65.Ta, 03.67.Ac, 03.65.Ud 



Nonlocality [l| and entanglement @] of composed sys- 
tems are distinguishing features of the quantum domain. 
Nevertheless, it is the possibility of locally broadcasting 
3 the state of a multiparticle system that broadly de- 
fines the nature of its correlations. Remarkably, even 
separable (nonentangled) states can be quantum corre- 
lated. This kind of quantumness has an important role 
that is not only related to fundamental physical aspects 
but also concerning applications in quantum information 
processing 043 and communication J3I Q , thermody- 
namics 1 91, quantum phase transitions [lfj, and biological 



systems [ll|. 



There are several unique aspects of quantum physics 
that discern it from classical theories. One of particular 
relevance to quantum information science is the impossi- 
bility of creating a perfect copy of an unknown quantum 
state 12[ . This fact is employed in some quantum crypto- 
graphic protocols [lH and, when extended to multipartite 
mixed states Q , can be used to classify the aspects of cor- 
relations in a composed system as classical or quantum. 
Let us consider a bipartite system described by the den- 
sity operator p and shared by parts a and b, with respec- 
tive Hilbert spaces H a and H b . The correlations in state 
p are said to be locally broadcast if there are auxiliary 
systems H ai ,rl a2 ,H bl ,H b2 and local operations (com- 
pletely positive, trace-preserving linear maps) A a : H a — > 
H ai ®ri a2 and A fc :rl h -> H bl (Z>H b2 such that p a ^ b ^ = 
A a ® A b (p) with I(p aibl ) = I(p a2b2 ), where the quantum 
mutual information— Z(p x «) = S{p x ) + S{p v ) - S(p xy ), 



Sip) 



-tr(plog 2 p), and p x — tr y (p xy ) — quantifies all 



the correlations between the systems x and y. It turns 
out that the correlations in a bipartite system can be lo- 
cally broadcast, and therefore have a classical nature, if 
and only if the system's state can be written as Q 
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where {|a,)} and are orthonormal basis for the 

subsystems state spaces W a and H b , respectively, and 
{Pi,j} is a probability distribution. 

The class of states in equation (Q]) is contained in the 
set of separable states — those states that can be gen- 
erated via local operations coordinated by communicat- 
ing classical bits — whose more general form is p sep = 
YliPiPi ® Pii where {p.;} is a probability distribution and 
Pi(Pi) is a valid density operator for the subsystem a(b). 
There are separable states that cannot be cast in terms 
of orthogonal local basis as those given in equation ([1]) 
and, therefore, present non-classical correlations. Un- 
derneath such states lies a non-classicality beyond the 
entanglement-separability paradigm, which can be quan- 
tified by a departure between classical and quantum ver- 
sions of information theory. One of the most popular 
measures for this kind of non-classicality is the quantum 
discord (l3l |. This quantifier has been receiving a great 
deal of attention EH- Arid it was proposed as a 
figure of merit for the quantum advantage in some com- 
putational models without or with little entanglement 

In general, measures of non-classical correlations in- 
volve complete knowledge of the system's state followed 
by extremization procedures. In the laboratory, the 
first task is implemented by quantum state tomographic 
methods and the second one is carried out by additional 
numerical manipulations. These procedures are demand- 
ing and propagate the unavoidable experimental errors. 
This observation motivates the search for alternatives 
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FIG. 1. (Color online) (a) Schematic representation of the op- 
eration sequence used to witness the non-classical nature of 
correlations, (b) Equivalent pulse sequence employed in our 
experiment. The thicker filled bars represent tt/2 pulses, the 
thinner bars indicate 7r/4 pulses, and the grey bars indicate 
7r/6 pulses with the phases as shown (negative pulse phases 
are described by a bar over the symbol). The pulses rep- 
resented as unfilled dashed bars are modified to achieve the 
different rotations necessary for the witness protocol. The 
dashed gradient pulse is applied to obtain the classically cor- 
related Bell-diagonal state. The time periods A, and jj 
represent free evolutions under the J coupling 



regarding the classification of correlations in quantum 
states. Once the nature of these correlations somehow de- 
termines what can and cannot be done with a given sys- 
tem, it is sometimes enough to know whether the correla- 
tions in that system have a classical or a quantum nature. 
To accomplish this last task it is convenient to have an 
observable witness for the quantumness of correlations in 
the system. However, as the state space of classical cor- 
related systems is not convex, a linear witness cannot be 
used in general, and we have to take advantage of a non- 
linear witness. For a wide class of two-qubit systems, p — 

(l ab + ELi ( A<r? ® l b + B l I a <g> a\ + Cial ® a\ )) /4, a 
sufficient condition for classicality of correlations is [l5[ 
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with O t = a? ® a\ for i = 1,2, 3 and 4 = Ei=i(^°i ® 
I fa + wd a ® erf). The af b) is the ith component of the 
Pauli operator in subsystem a(b). Ai,Bi, Z{,Wi G 3? 
with Zi , Wi randomly chosen and constrained such that 
Ei z f = Ei w i — 1- For the so-called Bell-diagonal class 
of states, p bd - (l ab + E-=i da? ® /4, W PM = is 



also a necessary condition for the absence of quantumness 
in the correlations of the composite system (in this case 



(Oa), 



[15(1). We can easily verify that the observ- 



ables in equation ([2]) can be written in terms of one com- 
ponent of the magnetization in one subsystem as (Oi) p — 
(a?®I%, with & = U a ^ h [R ni {e i )(>B) ni {9 i )]U a ^ 

where R ril (d l ) = R%.(9i) R^Wi), and R% h) (0i) is a 
local rotation by an angle Oi around direction n*, with 
Q\ = 0, 8-2 = 63 = 7r/2, n2 = y, and 713 = z. U a ^b is the 
controlled-NOT gate with the subsystem a as control. 
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FIG. 2. (Color online) (a) The 1 H spectra (normalized by 
the thermal equilibrium state spectrum) obtained after the 
witness circuit execution (with rotations R ni , Rn 2 , an d Rn 3 ), 
(b) witness expectation value, (c) quantum discord (light blue 
columns) and classical correlation (dark blue columns) mea- 
sured in three different initial states, pqc quantum correlated, 
pec classically correlated, and pr thermal equilibrium state. 
The dashed line represents the experimental error bound for 
determination of classically correlated (zero discord) states. 
The witness was measured directly performing the circuit de- 
picted in Fig I (a), while the classical correlation and the 
symmetric quantum discord was computed after full QST and 
numerical extremization procedures. The correlations are dis- 
played in units of (e 2 /ln2)bit [f8| . 



We experimentally implemented the aforementioned 
witness using the room temperature nuclear magnetic 
resonance (NMR) system. In this scenario the qubits 
(quantum bits) are encoded in nuclear spins and they 
are manipulated by radio-frequency (rf) pulses. Unitary 
operations are achieved by suitable choice of pulse am- 
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FIG. 3. (Color online) Real (left) and imaginary (right) parts 
of the deviation matrix elements reconstructed by QST for the 
two initial prepared Bell-diagonal states: (a) pqc quantum 
correlated (equivalent to Ci = 2e, C2 = 2e, and C3 = — 2e); (b) 
pec classically correlated (equivalent to C\ — 0, C2 = 0, and 
Cs = — 4e); and also for (c) the thermal equilibrium state. 
The deviation matrix elements are displayed in the usual com- 
putational basis, where |0) and |1) represent the eigenstates 
of a z for each qubit. The accuracy of prepared initial states 
can be estimated by the normalized trace distance from the 
ideal ones, 8(p idea i, p prep .)/e = tr|Ap Hea; — Ap prep . | /2 « 0.1 
(for both pqc and pec)- 



plitudes, phases and durations, and the transverse mag- 
netizations are obtained directly from the NMR signal 
[l6j |. The state of the two-qubit system is described 
by a density matrix in the high temperature expansion 
(where entanglement was ruled out), which takes the 
form p = I ah /4 + eAp, with e = hu L /4k B T ~ 1CT 5 
as the ratio between the magnetic and thermal energies 
and Ap as the deviation matrix 



12, UM. A carbon-13 



enriched chloroform (CHCI3) solution at 25° C was used 
in the experiments, with the two qubits being encoded 
in the 1 H and 13 C spin-1/2 nuclei. In order to experi- 
mentally demonstrate the witnessing protocol, two initial 
states were prepared by mapping them into the deviation 
matrix using the general pulse sequence scheme as shown 
in Fig. [TJ The first state corresponds to a quantum cor- 
related Bell-diagonal state, which is obtained from the 
thermal equilibrium by applying the pulse sequence for 
producing the pseudo-pure state |11), followed by the 



pulses that implement a pseudo-EPR gate [TjJ, see Fig. 
Q] [HI . The second state is a classically correlated Bell- 
diagonal state, obtained by applying a z-gradient pulse 



after the aforementioned pulse sequence. The witnessing 
of the thermal equilibrium state was also performed as a 
reference. The experimental procedure depicted in Fig. 
[T] was ran three times for each initial state in order to 
measure the magnetization (cr")^ in the states that 
leads to the two-point correlation functions (erf <g> <y\) p . 
So, the witness given in equation @ is directly measured 
(Fig.©. 
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FIG. 4. (Color online) Witness and correlations decoherence 
dynamics. The panel displays the measured witness and com- 
puted correlations for pqc [Fig. 3 (a)] relaxed during a time 
interval, t n = nSt (St = 55.7 ms, n = 0, 2, 11), before per- 
forming the witness measurement protocol. The red tick bars 
represent the witness expectation value (scale on the right), 
the dark grey section represents the amount of classical corre- 
lation, the light grey section represents the symmetric quan- 
tum discord, the entire grey bars (light and dark sections 
summed) display the quantum mutual information (scale on 
the left). The classicality bound is represented by the blue 
dotted line. The inset image shows the real part of the de- 
viation matrix elements reconstructed by QST for an inter- 
mediate classically correlated state. The effective transversal 
relaxation times are T 2 * = 0.31 s and T 2 * = 0.12 s, for X H 
and 13 C nuclei, respectively. The correlations are displayed 
in units of (e 2 /ln2)bit. 

For the sake of comparison, we performed a full quan- 
tum state tomography (QST) [19( of the initial states 
(displayed in Fig. [3]) and computed, from these data, the 
symmetric quantum discord p|, [2(| — Q(p ab ) — I{p ab ) — 
max| n „ n<f j_ I(x ab ), where l(x ab ) is the measurement- 
induced mutual information — and its classical counter- 
part [IH present in each state. The results are shown in 
Fig. [5J The correlation quantifiers are computed from 
the experimentally reconstructed deviation matrix in the 
leading order in e, following the approach introduced in 
Refs. [a, Since the error in the witness expectation 
value depends on many parameters (i.e., signal-to- noise 
ratio and residual rf pulse sequence imperfections), we 
used as a reference the thermal equilibrium state, which 
is supposed to have no correlations at room temperature 
6]. The witness measured for this state (W PT ) was about 
0.05, which is assumed to be the error margin for our ex- 
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periment. This introduces the bound shown in Fig. H](b) 
for a classically correlated (zero discord) state. 

The witness measured for the three initial states is dis- 
played in Fig. [5] (b). For the quantum correlated Bell- 
diagonal state the witness {W PQC ) is found to be about 
3.13 (far above the 0.05 bound), while for the classical 
correlated Bell-diagonal state (W pcc ) it is about 0.04, 
i.e., within the classicality cut-off limit. In fact, the wit- 
ness works perfectly in the present setup, in the sense 
that it easily sorts out quantum and classically correlated 
states. Figure [H (c) also displays the quantum discord 
computed from the experimentally reconstructed devia- 
tion matrices using the approach introduced in Ref. Q. 
As can be seen, the result for the quantum discord is in 
agreement with the witness, but the former is obtained 
after full QST and numerical extremization procedures. 
Finally, we followed the decoherence dynamics of the wit- 
ness, by letting the state pqc evolves freely during a time 
period t n , after this decoherent evolution we performed 
the witness circuit, and also a QST in order to com- 
pare the witness results with those for correlation quan- 
tifiers. The noise spin environment causes loss of phase 
relations among the energy eigenstates and exchange of 
energy between system and environment, resulting in re- 
laxation to a Gibbs ensemble. In Fig. 4 we observe, that, 
in the course of the witness and correlations evolution, 
the non-classicality is diminished until reaching an only 
classically correlated state. This occurs near the 1 H ef- 
fective transversal relaxation time. After such evolution 
period there are just reminiscent classical correlations, 
which arc also diminished resulting in an uncorrelated 
state (the room temperature thermal equilibrium state) 
after the spin-lattice relaxation time. Again, we obtain 
a fairly good agreement between the witness expectation 
values and the correlation quantifiers. 

Summarizing, we presented a direct experimental im- 
plementation of a witness for the quantumness of cor- 
relations (other than entanglement) in a composite sys- 
tem. Our work showed that it is possible to infer the na- 
ture of the correlations in a bipartite system performing 
only few local measurements over one of the subsystems 
(just three measurements for both pqc and pec)- The 
witness presented in Eq . (J5J) was generalized to higher- 
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dimensional systems [21[. Therefore, the methods em- 
ployed here can also be easily applied for witnessing cor- 
relations in systems with dimensions higher than two. 
Our strategy precludes the demanding tomographic state 
reconstruction and the numerical extremization methods 
included in quantum correlation quantifiers (like quan- 
tum discord). This method offers a versatile test-bed for 
the nature of a composite system that can be applied to 
other experimental physical contexts. Moreover, in such 
a proof of principle, we showed that non-classical correla- 
tions can be present even in highly mixed states as those 
in room temperature magnetic resonance experiments. 
The authors acknowledge financial support from 



Note added. — After the submission of this Letter, a re- 
lated work has appeared 22] , which employs other meth- 
ods to witness nonclassicality in an NMR system. 
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Supplemental Material - Details on experimental 
and calculation procedures 

NMR experiments 

Nuclear magnetic resonance (NMR) experiments were 
performed on a two-qubit system comprised by nuclear 
spins of 1 H and 13 C atoms in a carbon-13 enriched chlo- 
roform molecule (CHCI3). The sample was prepared by 
mixing 100 mg of 99 % 13 C-labelled CHC1 3 in 0.2 ml of 



where if [J-pj is the spin angular momentum operator in 
the a, j3 — x,y, z direction for 1 H ( 13 C); ip H (<f C ) defines 
the direction of the rf field and wf (wf) is the intensity 
of RF pulse for X H ( 13 C) nuclei. 

The first two terms describe the free precession of 1 H 
and 13 C nuclei about the external constant magnetic field 
Bo with frequencies u>h/2tt ks 500 MHz and ujc/^tt « 
125 MHz, respectively. The third term describes a scalar 
coupling between 1 H and 13 C nuclei at J w 215.1 Hz. 
The fourth and fifth terms represent the radio frequency 
(rf) field that may be applied to *H and 13 C , respectively. 
The 7r/2 pulse has a time setup of 7.4 /xs for 1 H and 9.6 
/is for 13 C at transmitter channel; or a time setup of 7.9 
/is for 1 H and 10.2 /xs for 13 C at decoupler channel. 

Spin-lattice relaxation times (Ti) of 2.5 s and 7 s for 1 H 
and 13 C, respectively, were measured using the inversion 
recovery pulse sequence. The recycle delay was set to 40 
s in all experiments. The effective transversal relaxation 
times, measured as the inverse of the spectrum line width 
at halt maximum, were found to be T 2 * = 0.31 s and 
T 2 * = 0.12 s for 1 H and 13 C nuclei, respectively. 

In order to experimentally demonstrate the witnessing 
protocol, two initial states were prepared by mapping 
them into the deviation density matrix using the general 
pulse sequence scheme shown in Fig. 1 of the main text. 
The first state corresponds to a quantum correlated Bell- 
diagonal state, which is obtained from the thermal equi- 
librium by applying the pulse sequence for producing the 
pseudo-pure state |11), followed by the pulses that imple- 
ment a pseudo EPR gate 0]. It is worth mentioning that, 
despite a true EPR gate would produce the same state, 
there is no need of using it in our propose, since the same 
kind of deviation density matrix can be produced with a 
pseudo EPR gate, with the advantage of using a smaller 
number of pulses. A classically correlated Bell-diagonal 
state, which is obtained by applying a z-gradient pulse 
after the aforementioned pulse sequence, was also pre- 



99.8 % CDCI3 in a 5 mm NMR tube. Both samples were 
provided by the Cambridge Isotope Laboratories - Inc. 
NMR experiments were performed at 25° C using a Var- 
ian 500 MHz Premium Shielded ( 1 H frequency), at the 
Brazilian Centre for Physics Research (CBPF). A Var- 
ian 5mm double resonance probe-head equipped with a 
magnetic field gradient coil was used. 

The rotating frame nuclear spin Hamiltonian account- 
ing for the relevant NMR interactions reads [l| 



(3) 

I 

pared. Finally, the witnessing of the thermal equilibrium 
state was performed as a reference for error estimation of 
the whole procedure. The y rotation in the witness proto- 
col was implemented directly by a single rf pulse, while 
for the z rotation the pulse sequence (§)_ -(§ ) a .-(f) y 
was applied in both nuclei. The controlled-NOT gate 
was achieved by the pulse sequence (f) C -U (^j)-(f) _ - 

(I )- y -(l )-,-(! )!;-(§ ) H y "(§ )" "(§ )^wlrere the super 
indices states for nucleus where the pulse is applied and 
U (jj) represents a free evolution under J coupling. In 
this pulse sequence that correctly implement a 

true CNOT gate, i.e., with the correct phases, are indeed 
necessary to produce the correct output for the witness. 
The protocol described in Rcfs. [4, 5] was used for quan- 
tum state tomography (QST). As described above, the 
experimental procedure depicted in Fig. 1 of the main 
text is ran three times for each initial state in order to 
measure the magnetization (<rf ^ in the states £j that 
leads to the two-point correlation functions (erf (g> 07 ) p . 
So, the witness given in equation (2) of the main text is 
directly measured. 



Measures of correlations 

The correlations in the high-mixed state find in the 
NMR context, p ab = P b /4 + sAp (e w 10~ 5 ), are com- 
puted by expanding the symmetric version of quantum 
discord @, @] 

Q(p ab ) = I(p ab ) - max l( X ab ), (4) 
{n ? ,m} 

where the quantum mutual information (expanded in the 
leading order in e), is given by 

I(p ab ) « ^ {2tr[(Ap ab ) 2 ] - tr [(Ap a ) 2 ] - tr [(Ap h ) 2 ] } , 



I 

H = -(uh- w#) If - (w c - If + 2tt Jlf If 

+wf (If cos tp H + If sin <p H )+uf (If cos ip c + 1% sin Lp c ) , 
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and the measurement-induced mutual information is tiers can be computed directly from the experimentally 

reconstructed deviation matrix and they were employed 

j, a b\ ^ _£ -f 2tr ffA ab ) 2 ] — tr fC A a ) 2 l — tr [(A b ) 2 H to observe the quantumness of a room temperature NMR 

^ X J ~ In 2 * r Ll X ^ J r U X J J r H X J J J- , quadrupolar S y S tem 0. 



with x ab = I ab /4 + eAx ab as the state obtained from 
pab through a complete projective measurement map 

(A x ab = Eij-n? ® n b (Ap ab )nf g) n b ). Ap Q ( b ) = 

tr b ( a ){Ap ab } is the reduced deviation matrix while Ax°^ b ' 
stands for the reduced measured deviation matrix in the 
subspace a(b). The classical counterpart of Eq. 0] is 
C{p AB ) — max| na n6 | l(x AB )- It is worth mentioning 

that Q(p ab ) — if and only if p ab can be cast in terms 
of local orthogonal basis. In other words, the symmetric 
quantum discord is zero if and only if p ab has only clas- 
sical correlations or no correlations at all (in this case 
the correlations present in p ab can be locally broadcast 
Q). The aforementioned symmetric correlation quanti- 
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